Perturbative predictions for B c meson production in 

hadronic collisions 



Marco Masetti 

Dipartimento di Fisica, Universita di Roma "Tor Vergata", 
and I.N.F.N., Sezione di Roma II, Viale della Ricerca Scientiflca 
1-00133 Roma, Italy 

Francesca Sartogo 

Dipartimento di Fisica, Universita di Roma "la Sapienza", 
and I.N.F.N., Sezione di Roma, Piazzale A. Moro 2, 
1-00185 Roma, Italy 

March 31, 1995 

Abstract 

Perturbative cross section for direct B c meson production in gluon gluon scatter- 
ing gg B+bc is calculated and compared with other existing results. Predictions 
for hadronic B c production at Tevatron and LHC are presented and the main sources 
of uncertainties are discussed. 
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In the last years many theoretical studies on B c meson properties were published, but 
an experimental identification of this pseudoscalar bound state formed by two quarks with 
different flavours is still missing, due to the very small production cross section. 

The mass spectrum of bc bound states is expected to be similar to the spectra of other 
mesons with two heavy quarks (charmonium and bottomonium) and predictions based on 
potential models and on QCD sum rules [1, 2] indicate the presence of a rich structure of 
excited states below the fragmentation threshold in a B and a D meson. There is however 
an important difference with respect to unflavoured quarkonia: the excited states cannot 
decay directly into light mesons, but decay only to the pseudoscalar B c ground state, 
due to the flavour conservation in strong and electromagnetic interactions. The lifetime 
of the weak-decaying B c meson is expected to be in the range of D mesons lifetimes 
[3, 4]. Exclusive B c decays were studied by several authors [3, 4, 5]. Some decay modes 
are interesting for experimental detection, for example the semileptonic decay with three 
charged leptons in the final state (B c — > J/ip + I + v\ followed by the electromagnetic 
decay J/ip — > l l+ + with composite branching ratio of about 10~ 2 , and the decay in 
J/ip + 7r, with branching ratio of the order of 10~ 3 . The lack of experimental detection of 
B c mesons is a consequence of the small production cross section: in fact the production 
via electromagnetic or strong interaction requires two additional heavy quarks in the final 
state, and is therefore suppressed at typical experimental energies. The production of a 
B c without other heavy particles should be possible in processes like e + u — > W* — > B c , 
which however involve weak interactions and are strongly suppressed due to the high 
virtuality of the W*. 

The standard estimates of B c production cross section are based on perturbative cal- 
culations with the heavy quark bound state described in the zero binding energy limit. 
The simplest but important case of B c production in e + e~ annihilation involves only four 
Feynman diagrams and was studied in detail [6]. The total branching ratio for Z° with a 
B c + or a B c ~ in the final state (including the contribution due to the production of excited 
B c states which subsequently decay to the ground state) is predicted to be 140 • 10~ 6 i. e. 
a(Z° -> X) ~ 10-V(Z° -> b b X). An estimate obt ained with the Montecarlo event 
generator Herwig, which for Z° decays in B c ^ gives similar results, indicates that also at 
high energy hadronic colliders like Tevatron and LHC the B c ^ production cross section 
would be around 10~ 3 times the bb production cross section [7]. 

In the perturbative approach the relevant process for B c hadroproduction is the gluon 
gluon scattering g g — >■ B+ b c, which at the lowest order involves 36 diagrams. The 
process qq — > B+ b c is expected to give a negligible contribution to the hadronic B c 
production cross section. The calculation is extremely simplified if one considers only the 
fragmentation terms [8], which were expected to dominate at high energy; however this 
approach fails at least in the photoproduction process, as shown in [9]. The complete 
perturbative calculation for g g — >■ B+ b c was performed by several authors [11, 12, 13], 
with results in disagreement with each other. 

From this summary of theoretical predictions about the B c mesons it should be clear 
that high energy hadronic colliders give the best chances to identify these particles, and it 
is therefore necessary to reduce the large uncertainties due to the discrepancies between 
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existing predictions. For this reason in this letter we present the results of an independent 
study on direct B c hadroproduction. The dependence of the results on different choices 
of the gluon distribution functions and on the scale for a s , which are the main sources 
of uncertainty, is discussed. Our results are definitely incompatible with [11] and [13], 
which however give predictions differing with each other at least by a factor 20, difficult 
to ascribe only to different choices of gluon distribution functions. On the other hand our 
results are essentially in agreement with [12]. 

The calculation technique is standard: in the approximation of negligible binding 
energy and relative momentum in the bound state, which is justified for bound states with 
only heavy constituents, the quarks are on shell and have parallel momenta (m^p^ = m c p£) 
and the meson mass is Mb c = m c + ra&. The 36 Feynman diagrams for the process 
g g — >■ B+ b c can be obtained from the 13 diagrams of fig. 1 performing all possible 
interchanges of initial gluon momenta and of final quark flavours. The rule to describe 
the bound state in this approximation is to insert as 'bound state vertex' the spin projector 

S cc ' fn 

TB < T(jf> + M Bc ) 



where c and c' are colour indices and T = 75 for a pseudoscalar meson while r = jk for a 
vector meson. The decay constant fs c is defined by 

< 0\h^ 5 c\B c (P) >= if Bc P^ 

and in the non relativistic limit it is proportional to the wave function in the origin: 



>«■ = ^'* (0)| - 

The simple substitution of the sum on external gluon polarizations X)»=i 2 £ iii £ iv with g^ v 
would give ghost contributions. It is possible to avoid ghosts summing only on physical 
transversal polarizations: 

Ek^Tt^ ~\~ k-^Tl^ k^kjj 
£ ifj. £ iu 9l*v T 



= 1,2 



k ■ n {k ■ to) 2 



where k is the gluon momentum and we choose to m = {k a ^ + k^/^fi. The number of 
independent terms in the cross section can be greatly reduced using relations between 
amplitudes connected by exchange of the initial gluons or of final quarks flavours and 
momenta. Using the algebraic program Schoonschip [14] we calculate the squared ampli- 
tude and integrate analytically to obtain the differential partonic cross section d 2 a / ' dpxdy, 
where px is the transverse momentum and y is the cm. rapidity of the produced B c . The 
subsequent integrations are performed both using the Montecarlo integration program Ve- 
gas [15] and using simple algorithms for numerical evaluation of low dimensional integral, 
to check the convergence of the result. 
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VI(GeV) 


a 


«x[ll] 


a [12] 


20 


23.6 


5.3 


22.5 


30 


26.0 


9.2 


— 


40 


22.3 




30.8 


60 


15.5 


8.6 




80 


11.2 






100 


8.4 




16.2 



Table 1: Cross section (in pb) for g g — >■ B+ b c. All the results are rescaled to f Bc = 500 
MeV. 



As a preliminary step and a check of our calculation we studied the direct production 
in photon photon collisions 77-* b c. At the lowest order of perturbation theory this 
process involves 20 diagrams which can be obtained from those of fig. 1 eliminating the 
diagrams containing non abelian vertices, and replacing the external gluons with photons. 
This calculation was already performed in [9] and [10]. Our results for this process, using 
the same values for the parameters ra&, ra c , f Bc , a s , a em , are in good agreement with 
those of [9] and [10]. 

For the process g g — >■ b c we use the following values of the parameters: 

f Bc = 500 MeV, a s = 0.2, M Bc = 6.3 GeV, m b /m c = 3 

Note that the partonic cross section is a function of the kinematical variables times f Bc a s 4 , 
and therefore the rescaling of the results for different values of f Bc and a s is easy to obtain. 

The total cross section for this process as a function of the energy in the partonic center 
of mass is plotted in fig. 2, while the differential partonic cross sections da/ dpT and 
da/ dy are shown for different values of s in fig. 3 and fig. 4. 

In Tab. 1 the total partonic cross section for several values of s is compared with the 
analogous existing calculations [11, 12], rescaling all the results to the same value of f Bc . 
Unfortunately the comparison is made difficult by the small amount of information on 
the partonic cross section before the convolution with the gluon distribution functions. 
The results of [11] are significantly smaller than ours. Even if it is not clear to us if their 
results refer to a fixed value of a s or a running one, and in this last case which scale they 
choose, the difference cannot be explained simply by the running of a s . 

In [12] a s is fixed to 0.2 and f Bc = 570 MeV; with the appropriate rescaling of f Bc our 
results are in reasonable agreement with theirs in the low s region. Given the behaviour 
of the partonic cross section, the discrepancies between the values we obtain and those 
obtained by [12] at higher energies are not really relevant in the evaluation of the hadronic 
cross section, after the convolution with the gluon distribution functions. Finally, a com- 
parison with [13] is impossible at this level, since only results on the hadronic cross section 
are discussed there. 
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gluon 
distribution 


scale 


a (nb) 


a (nb) 
this paper 


0.12 
1.8 
16. 


[11] 


EHLQ 


a/4 


4.5 10" 3 
1.0 
12.3 


60. 10" 3 
8.5 

78 


1.8 
16. 


[12] 


[12] 


8 * 
8 * 


9.4 

151 


8.8 
140 


1.8 
1.8 
1.8 


[13] 


CTEQ 


S 

a/4 
4M Bc 2 


18.3 
29.6 
31.5 


2.5 
3.8 
3.9 



Table 2: Comparison of hadronic cross section (in nb). All results are rescaled to fs c = 
500 MeV. In the fourth column we indicate the scale for gluon distribution functions and 
running a s ; the * indicates results obtained for fixed a s = 0.2, evolving only the gluon 
distribution. In the fifth and sixth column the cross sections obtained in the cited papers 
and the values we obtain using the same prescriptions are reported. 



The cross section for the hadronic processes p p — > B+X and p p — > B+X can be 
obtained by the convolution of the partonic cross section with the gluon distribution 
functions g(x,Q 2 ): 

cr(a) = ( dxi f dx 2 g(x 1 , Q 2 )g(x 2 , Q 2 )9(x 1 x 2 s - AM Bc 2 )cr{s = x 1 x 2 s). 
Jo Jo 

In figs. 5 and 6 the differential cross section in px and yi a b obtained with the MRS(A) 
gluon distribution function [16] is plotted for the energies corresponding to Tevatron (1.8 
TeV) and LHC (14 TeV). 

In Tab. 2 we compare the values for direct B c + production cross section in proton pro- 
ton collisions reported in [11, 12, 13], with the results we obtain using the same prescrip- 
tions and gluon distribution functions (EHLQ [18], CTEQ [19] and the parametrization 
described in [12]). In [11] the EHLQ gluon distribution function is used fixing the evolu- 
tion scale Q 2 at J/4. With the same gluon distribution and the same Q 2 we obtain for 
the hadronic cross section values greater of those of [11] especially at the lower energies. 
This comparison confirms the disagreement with [11] already observed for the partonic 
cross section. 

To compare our predictions with [12], we fix a s = 0.2 as they do, and we use the same 
parametrization of gluon distribution function they describe in [12] with an evolution 
scale Q 2 = s. There is agreement between our results and those of [12], especially at the 
lower energy. The results of [13] are much larger than all other existing predictions. We 
obtain results which are smaller than those of [13] by at least a factor 7, using the same 
gluon distribution (CTEQ) and the same evolution scale as they do, and therefore we 
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3.8 
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5.8 


7.8 
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84 
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9.4 


60 
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Table 3: Hadronic cross section (in nb) with f Bc = 500 MeV; Q 2 is the scale for gluon 
distribution functions; a and a are the cross sections calculated with running a s (Q 2 ) and 
with fixed a s = 0.2 respectively. We put in evidence the maximum and minimum value 
obtained for running ct s (Q 2 ) using the more updated gluon distribution parametrizations. 



conclude that the incompatibility of our results with those of [13] is not related to the 
choice of gluon distribution functions and can be only a consequence of a disagreement in 
the calculation of the partonic cross section a(g g — >■ b c). 

In Tab. 3 the results we obtain with several gluon distribution functions (MRS(A) and 
MRS(G) [16], GRV [17], EHLQ [18], and the parametrization of [12]) are shown. While 
the choice of the gluon distribution function parametrization EHLQ and [12] is mainly 
dictated by the comparison with other existing results, it is also interesting to show the 
values obtained with other updated sets of gluon distribution functions. In order to discuss 
the possible dependence of the results on the evolution scales for the gluon distributions 
and for the running coupling constant a s , we vary Q 2 between \M B 2 and s. The case of 
fixed a s (a s = 0.2) is also considered, as suggested in [12]. We will not attempt to justify 
a particular choice of the evolution scale in the leading order calculation described here; 
it is however clear that results obtained fixing a s = 0.2 must be considered as optimistic 
higher limits on the predicted cross sections. It should be noted that the relatively high 
result one obtains using EHLQ and a s running is due to the relatively high value of Kqcd 
rather than to the shape of the gluon distribution function. The difference in the choice of 
the gluon distribution function however emerges more clearly at the energy corresponding 
to LHC. While at the energy of Tevatron the relevant x region is x > 10~ 4 , in the case 
of LHC it becomes x ^ 10~ 6 . All gluon distribution parametrizations are expected to be 
valid only above a given x m i n (a typical value is x m i n = 10~ 5 ) and therefore the cross 
sections at higher energies are typically underestimated. The only exception is the GRV 
gluon distribution function, for which x m i n = 10~ 6 , but it has to be noted that the shape 
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of the GRV gluon is considerably steeper than, for example, the MRS(G) gluon, which 
comes from a recent analysis of the latest measurements of the structure function F 2 at 
HERA [16]. The predicted cross section for direct B c production at Tevatron and LHC, 
taking into account the uncertainties on the knowledge of the behaviour of the gluon at 
very small a;, can be obtained with MRS and GRV gluon distributions: 

vxrectipp -+ BtX) = {m f ^ eV)2 (3-2 - 6.4) nb (Tevatron; 1800 GeV); 



<7 d irect( P P -+ BtX) = {m f ^ eV)2 (0-04 - 0.12) /xb (LHC; 14 TeV). 

A possible important source of uncertainty in the predicted cross sections is the fac- 
tor /| c . The existing predictions vary from f Bc = 160 MeV [20] to f Bc = 600 MeV 
[21]. Varying f Bc between so different values would give a contribution to the uncertainty 
much larger than those due to the choices of parton distributions and of evolution scales. 
A more careful inspection of the predictions shows that potential models have a much 
smaller uncertainty {fs c = 500 ± 80MeV: for a detailed discussion and review of predic- 
tions for f Bc see [2]), while QCD sum rules give results spread in a wide interval, due to 
the intrinsic ambiguities of the method. However in [2] it is asserted that the most reliable 
predictions from QCD sum rules are compatible with potential models {fs c = 460 ± 60 
MeV). The simpler way to obtain predictions for f Bc and for the cb spectrum is to rescale 
the charmonium and bottomonium spectra to the intermediate case of cb bound states: 
this is the role of potential models which can be considered as a sort of parametrization 
of quarkonium spectra which gives predictions for B c . For this reason it seems to us that 
potential models are presently more reliable than QCD sum rules for phenomenological 
predictions on cb spectrum and wave functions. Considering only potential model predic- 
tions, the uncertainty in B c production cross sections due to the factor f B is about 30%, 
much smaller than the uncertainty from a s and gluon distributions. 

Finally we want to give an estimate of the total B c production cross section in hadronic 
collisions, including the contribution due to the production of excited cb states which sub- 
sequently decay to pseudoscalar B c . The main contribution should come from vector B* 
production. In [12] this cross section is calculated obtaining <r(pp — > B*X) ~ 2.5a(pp — > 
B C X). Neglecting further contributions from other excited B c states we conclude that the 
total cross section for the production of B^ and B~ would be 

<Ttotai(pp ->■ B+X) ~ 20 - 50 nb (Tevatron; 1800 GeV); 

<rtotai(pp -> B+X) ~ 0.3 - 0.8 /xb (LHC; 14 TeV). 
This correspond to at least one thousandth of the total bb production cross section. 
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In this letter we discussed the main sources of uncertainties in the calculation of the 
cross section for direct B c meson production in hadronic collisions at the lowest order 
of perturbation theory and in the zero binding energy limit. The uncertainties related 
to the choice of input parameters and of gluon distribution parametrizations amount 
roughly to a factor 3, too small to explain the large discrepancies between the different 
existing calculations. Finally we gave an estimate of the total B c production cross section, 
including the contribution of excited states, for Tevatron and LHC. 

We thank M. Lusignoli for having drawn our attention to the problem treated in this 
letter and for many useful discussions and comments. We are also grateful to A. K. 
Likhoded, V. Lubicz, S. Petrarca and B. Taglienti. 
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Figure Captions 



[1] Feynman diagrams for g g — >■ B+ b c. The complete set of 36 diagram can be 
obtained from the diagrams shown in the figure performing all possible interchanges 
of initial gluon momenta and of final quark flavours. 

[2] Total partonic cross section a(g g — > B+ b c) (in pb). 

[3] Partonic p T distribution da/dp T (in pb GeV" 1 ) for ^/l = 20, 30, 40, 60, 80 and 100 
GeV (in order of increasing px endpoint). 

[4] Partonic rapidity distribution da/dy (in pb) for = 20, 30, 40, 60, 80 and 100 
GeV (in order of increasing rapidity endpoint). 

[5] Differential cross section da/dpT (in nb GeV -1 ) for direct B c production calculated 
with MRS(A) gluon distribution function for LHC {^/s = 14 TeV) with the evolution 
scale Q 2 = 4:M^ c (upper solid line) and Q 2 = s (lower solid line), and for Tevatron 
(y^ = 1.8 TeV) with the evolution scale Q 2 = 4:Mg c (upper dotted line) and Q 2 = s 
(lower dotted line). 

[6] Differential cross section da/dyi a b (in nb) for direct B c production calculated with 
MRS(A) gluon distribution function for LHC (y^ = 14 TeV) with the evolution 
scale Q 2 = 4:M^ c (upper solid line) and Q 2 = s (lower solid line), and for Tevatron 
(y^I = 1.8 TeV) with the evolution scale Q 2 = 4:M^ c (upper dotted line) and Q 2 = s 
(lower dotted line). 
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Figure 1: 
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Figure 5: 
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Figure 6: 
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